Influence of classical resonances on chaotic tunnelling 
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Dynamical tunnelling between symmetry-related stable modes is studied in the periodically driven 
pendulum. We present strong evidence that the tunnelling process is governed by nonlinear reso- 
nances that manifest within the regular phase-space islands on which the stable modes are localized. 
By means of a quantitative numerical study of the corresponding Floquet problem, we identify the 
trace of such resonances not only in the level splittings between near-degenerate quantum states, 
where they lead to prominent plateau structures, but also in overlap matrix elements of the Floquet 
eigenstates, which reveal characteristic sequences of avoided crossings in the Floquet spectrum. 
The semiclassical theory of resonance-assisted tunnelling yields good overall agreement with the 
quantum-tunnelling rates, and indicates that partial barriers within the chaos might play a promi- 
nent role. 
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I. INTRODUCTION 

Since the very first quantitative studies of tunnelling 
in a chaotic system 0, > it has been clear that a minute 
scrutiny of the associated classical dynamics was required 
in order to understand even the most coarse features of 
the quantum behaviour of such systems. In essence, tun- 
nelling is a semiclassical concept since it refers to a quan- 
tum process — typically a decay or the oscillation of an 
averaged observable — that is forbidden at a classical 
level. But, despite the numerous successes of semiclas- 
sical computations in quantum chaos, the questions of 
which and how classical objects can be used to under- 
stand tunnelling and to compute, say, their characteris- 
tic time scales, have been remaining widely open for six- 
teen years. The stakes in the battle are important since 
non-integrability is the generic rule of multidimensional 
systems and tunnelling may play a crucial role in their 
transport properties. Moreover, it is one of its signatures 
that chaotic tunnelling can be modified on several orders 
of magnitude by the slightest variation of any classical 
or quantum parameter; therefore a deep understanding 
of chaotic tunnelling is required to control the process, 
what may be an advantage in delicate quantum experi- 
ments and, hopefully, give rise to an extremely sensitive 
quantum tool. Some promising clues have been provided 
in this direction by numerical studies and experiments 
with cold atoms ^ |B< i^J but also with microwave cav- 
ities 0, where tunnelling signatures for processes that 
are forbidden by ray optics were observed. 

Of course, in order to capture the typical exponentially 
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small tunnelling effects, it is expected that the classical 
dynamics should be complexified. It is well-known ,.9, .loj 
that the complex solutions of Hamilton's equations are 
actually involved in the interpretation of tunnelling of 
autonomous systems with one degree of freedom. In the 
mid 90's, the first observations [Hi 0113, 113 that com- 
plex periodic orbits allow one to reproduce quantitatively 
some feature of chaotic tunnelling gave hope that a semi- 
classical strategy was indeed possible, even though the 
complexified classical tori are generically destroyed in 
chaotic systems But to deal with tractable semi- 

classical trace formulae a la Gutzwiller, a general crite- 
rion for selecting the complex periodic orbits was still 
lacking; this need became an emergency when it was un- 
expectedly discovered [l6l IT^ that chaos reveals itself in 
the complex phase-space through some fractal structures, 
the so-called Laputa islands, that look like agglomerates 
of complex trajectories. It is only recently that some 
encouraging significant steps were done for retaining the 
relevant semiclassical skeleton . A lot of work remains 
to be done in that direction, especially if one wants to 
deal with continuous systems where time can (and must) 
be complexified as well, unlike what occurs in discrete 
maps. 

The second strategy to cope with chaotic tunnelling is 
not purely semiclassical but rather calls up random ma- 
trix theory. Since the seminal work presented in Ref . |l9| , 
a fruitful approach of quantum chaos is to replace a 
chaotic but deterministic Hamiltonian by a random ele- 
ment of an ensemble of matrices that only encapsulates 
the global symmetries. These hybrid techniques, with 
both semiclassical and statistical ingredients, first al- 
lowed us to qualitatively understand the so-called chaos- 
assisted tunnelling, i.e the observation 20] that tun- 
nelling is increased on average as the transport through 
chaotic regions grows [2l[|22,|23. However, the extreme 
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sensitivity of tunnelling renders the predictions very dif- 
ficult even if just an order of magnitude is required. 
This is also true for the seemingly simple case of near- 
integrable dynamics where it was shown, on a discrete 
quasi-integrable quantum map, that the internal reso- 
nances may enhance the transitions by several orders of 
magnitudes ^2l, . This resonance-assisted tunnelling 
is also at work in_discrete systems where chaos is much 
more developed 0. The aim of the present paper is to 
show that the ideas in Refs. [2^ H^l are strengthened, 
now in a continuous system, by a systematic analysis of 
level dynamics and the phase-space representation of the 
quantum states. 

We shall begin in section^ with a short presentation 
of the general framework of chaotic tunnelling and the 
model we choose in order to study it. We will be con- 
cerned with a typical signature of tunnelling, namely the 
period of Rabi oscillations between two wells that are 
separated by a dynamical barrier. Then we will show 
(section IIII|) with a simple argument, that the attempt 
to re pro duce the average tunnelling periods presented in 
Ref. |23| is far from being complete precisely because it 
ignores the resonances, among other things. In section llVl 
we will give an illuminating illustration of a characteris- 
tic feature of chaotic tunnelling : it appears to be 
a collective effect in level dynamics where not just one 
third state crosses a tunnelling doublet. We will give a 
phase-space picture of resonance-assisted tunnelling and 
confirm, that taking into account the resonances is un- 
avoidable if we want to reproduce or predict the ave- 
rage behaviour of tunnelling transitions. In section 
we show that the ideas of Refs. [2^ actually provide 
a good estimate for the average tunnelling rate in our 
model. 



II. GENERAL FRAMEWORK OF CHAOTIC 
TUNNELLING 

The simplest non-integrable Hamiltonian models are 
either time-dependent one-dimensional (ID) systems or, 
equivalently, autonomous systems with two degrees of 
freedom where the Hamiltonian is the only constant 
of motion. Seen from the classical point of view, a 
generic potential induces a cascade of non-linear resonan- 
ces whose overlap generates chaos [s^, |^ . One minimal 
continuous model that encapsulates these typical prop- 
erties is a ID time-dependent system whose Hamiltonian 
is 
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H{p,q;t) = — cos{q + 2Trt/T)-— cos{q-2-Kt/T) , 

(1) 

where p and q denote canonical action-angle variables 
that are coordinates on a phase-space having the topol- 
ogy of a cylinder: it is 27r-periodic in position (angle) q 
and infinitely extended in the momentum (action) direc- 
tion. The model ^ can be seen as the most natural 
normal form where we keep only two overlapping reso- 



nances i.e. the first two time-dependent Fourier compo- 
nents of the potential. The Hamiltonian H can also be 
interpreted, provided a change of frame is performed, as 
a pendulum driven by a periodic wave. Moreover, it is 
actually very similar to the effective Hamiltonian that 
can be realised in experiments on cold atoms 32]. We 
shall consider the period r = 27r in the following. 

The two parameters 7^ control the size of the two sta- 
ble islands I± located in phase-space near p = ±1 respec- 
tively. When increasing the 7's from zero, we leave the 
quasi-integrable free-like motion and rapidly (at 7 0.2) 
reach a mixed regime where the two stable resonance is- 
lands are fully surrounded by a chaotic sea. If an initial 
condition lies inside one island, the classical motion will 
remain trapped forever within a very thin quasi-one di- 
mensional layer without the possibility of escape (this 
evasion is the forbidden process that tunnelling will be 
concerned with ^S^). Alternatively, a trajectory starting 
in between the two islands is chaotic: without being fully 
ergodic since it cannot enter the islands, the absence of 
a sufficiently large number of constants of motion allows 
the system to explore large areas and to develop an ex- 
ponential sensitivity on initial conditions. Still increas- 
ing the 7's, the cascade of secondary resonances inside 
the islands erode gradually the regular zones and make 
them dissolve completely in the chaotic sea. Figure Q] 
shows some Poincare sections of the classical dynamics, 
i. e views taken stroboscopically at every integer multiple 
of the time period t = 27r of (^). 

At a quantum level, the quantization of a Hamil- 
tonian like has, for a long time, been employed 
as a very natural continuous model for studying tun- 
nelling |34li,35. 36,J. The time-dependence is implemented 
within the Floquet theory where the phases e of the eigen- 
values of the evolution operator at time r play the role of 
the energies for autonomous systems. The eigenvectors 
associated with these so-called quasi-energies e will be 
called the (eigen)states of the system. When h is small 
enough compared to the typical size of the stable islands, 
some states, the regular states {|</'n)}ri6{o,i,... }j appear to 
have their Husimi distribution localised on the two main 
stable islands while other states, the chaotic states, are 
found to be delocalised in the chaotic sea. 

Following the same route that is explained in great 
detail in Ref. [2^ §§ II & HI], we will keep a discrete 
two-fold symmetry, namely the time-reversal symmetry, 
by considering 7+ = 7_ = 7 throughout this paper [33. 
This allows us to clearly identify the tunnelling process 
through the existence of small energy scales, the dou- 
blet splittings Ae„ , which are associated with large time 
scales, namely the periods 27rfi,/Ae„ 3> t of tunnelling 
oscillations between the stable islands. In other words, 
the magnitude of the splittings Ae„ = Je^ — e^J between 
the states that are symmetric {\4>n)) and antisymmetric 
(](/)")) with respect to the time-reversal symmetry mea- 
sure the importance of tunnelling between the two stable 
islands. The actual challenge of chaotic tunnelling is to 
understand and hopefully predict the behaviour of Ae„ as 
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FIG. 1: Some Poincare sections (r = 27r-stroboscopic plots) 
corresponding to Hamiltonian Q for several values of 7+ = 
7_ = 7. Increasing 7, between 7 ~ 0.1 and 7 ~ 0.3, the 
two resonances at (p, (?) — (±1,0) start to significantly over- 
lap, and the systems enter in a mixed regime where both 
regular and chaotic dynamics coexist. From 7 ~ 0.3 to 
7 ^ 0.8, the two stable islands T± are completely surrounded 
by a chaotic sea, and although they are related by the time- 
reversal symmetry, they are disconnected from each other at 
a classical level, since dynamical barriers prevent the real 
classical trajectories starting in one island to escape from 
it. For 7 ~ 0.72, T± have developed a relatively wide 3/7 
resonance, which is not the case for 7 ~ 0.67 where the cor- 
responding 3/7 resonance chain is hardly visible. The insets 
in FIG.|21show magnifications of the islands in these last two 
cases. 



a classical (7) or quantum (h) parameter is varied within 
a regime where the quantum scales are small enough to 
resolve the classical scales of the chaotic structures. 

Unlike what occurs for integrable multidimensional 
systems where Ae„ is a smooth monotonic function given 
by Ae oc |3^ with A being a typical classical ac- 

tion that characterizes the tunnelling barriers, the split- 



tings can display huge fluctuations in the mixed regular- 
chaotic case, which were soon identified as a signature 
of chaos [22, 13 ■ -But as far as isolated fluctuations are 
observed this point of view must be amended. One iso- 
lated fluctuation is clearly associated with a third state 
whose (quasi)energy nearly degenerates with the dou- 
blet [SJ, I33 as 7 or ?i is varied. Chaos is not necessa- 
rily relevant here since such fluctuations can be observed 
when the chaotic layers are too small (and h too large) to 
be resolved by the quantum waves. It may happen that 
the third state is also a regular state localised on another 
regular EBK torus, possibly belonging to another stable 
island. What deserves the name of chaotic tunnelling is 
a radical change of regime where the fluctuations are not 
isolated anymore and where the coarse-grained behaviour 
of Ae does not follow a monotonic law. It was shown in 
Ref. § V] , that there actually exists a rather abrupt 
transition between a quasi-integrable tunnelling regime 
and a chaotic regime and that this transition occurs pre- 
cisely when the quantum eyes can resolve the chaotic 
classical structures. 

While the study of statistics of the splittin gs i s neces- 
sary to give an insight into chaotic tunnelling j23H39l|. we 
will focus here on the average behaviour of the splittings. 
It is far from obvious that we can conceptually justify 
the distinction between "large" scales (the average) and 
"small" scales (the fluctuations) in the variations of Ae. 
Most probably, there may exist a whole hierarchy of such 
variations. From the numerical point of view, however, 
there seems to be an overall modulation of Ae, both ob- 
served in maps (|SfIG. 3], [H FIGs. 1 & 2]) and in 
continuous systems p, FIG. 7]. This modulation is pre- 
cisely the very object of Ref. |26| (see also Ref. |27j|) and 
of the present paper. 



III. KEY ROLE OF RESONANCES 

In the general context recalled in the previous section, 
where tunnelling transitions occur between two symme- 
tric but disconnected stable classical islands, the follow- 
ing estimate has been proposed for a typical chaotic tun- 
nelling splitting [2^ eq. (4)]: 



Ae ~ Ml 



r(27V, 4iV) 
"r(2iV + 1;0) 



(2) 



Here, T stands for the incomplete Gamma function, 
N = A/{2nh) denotes the semiclassical estimate of the 
number of states localised in one island of area A, and 
il represents an unknown prefactor, with the dimension 
of an inverse time scale, which does not depend on h. 
Though the origin of the formula remains obscure as 
Eq. (O is not explicitly proven by their authors, it gave, 
in Ref. [2^, good agreement with numerical computa- 
tions provided we are ready to accept an unreasonably 
large ambiguity on the unspecified proportionality fac- 
tor: for ~ 40, the estimate ^ varies by five orders 
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of magnitude (10""'^ versus 10~^) in the two cases con- 
sidered in FIG. 2 of Ref. where the areas A of the 
stable islands are of the same order. More recent calcu- 
lations within the "kicked Harper" model, however, re- 
vealed substantial deviations between Eq. (0) and the ex- 
act quantum tunnelling rates in the semiclassical regime, 
where a reasonable agreement was only found in the deep 
quantum limit of large h . 




FIG. 2: (Color online) Comparison of the tunnelling split- 
tings Aeo for 7 = 0.72 (black line) and 7 — 0.67 (gray/red 
line). The corresponding stable islands, magnified in both 
cases with the same scaling factor, are shown in the insets 
(see FIG.Qc) and d) for a full Poincare surface of section). 
The exponentially large difference between the two curves is 
explained by the different internal structure of the islands: 
for 7 = 0.72, the 3/7 resonance chain is much more developed 
than for 7 = 0.67. 

Indeed, no good estimate of Ae can be obtained if the 
only classical parameter on which the theory depends is 
the area A of the island. As shown in FIG. 13 the in- 
ternal classical structure of the islands must be involved 
in one way or another. We have plotted here the split- 
ting Aeo between the two "central" states \4>^) localised 
in the two islands; more precisely, they both were se- 
lected by the criterion of having the maximal overlap 
with a coherent state that is located on the central sta- 
ble periodic orbit of period t. The two graphs. Acq as 
a function of 1/h, are shown for the two classical pa- 
rameters 7 — 0.67 and 7 = 0.72 (see also FIGQ c) 
and d)). In these cases, the stable symmetric islands T± 
have an area A of the same order of magnitude, but ex- 
hibit a rather distinct internal structure: for 7 = 0.72, 
T± have developed a wide 3/7-resonance chain compared 



to the 7 = 0.67 case. In the log-plot, the two graphs 
strongly differ. Not only the fluctuations of Acq hardly 
match, but also the average behaviour is completely dis- 
tinct in a semiclassical regime where > 10. A dis- 
crepancy of about five orders of magnitude can be clearly 
observed for ~ 16: for 7 — 0.72, we have Aeg ^ 10^^ 
compared to Acq ~ 10"" for 7 0.67. At 1/ft - 25, 
on the other hand, the splittings for 7 = 0.72 are about 
10"'^ smaller than for 7 = 0.67. 

In any case, no estimate formula Aeo{A,h), which 
would produce approximately the same graphs for 7 = 
0.67 and 7 — 0.72, can be satisfactory. At this stage 
it is more than plausible that the internal structure 
of X± must be taken into account in any attempt to es- 
timate Ae„ with the help of classical ingredients. 



IV. MULTILEVEL CROSSING 

To obtain more insight, we analyze in this section the 
"level dynamics", i.e. the changes of the quasi-energy 
spectrum as 1/fi. or 7 is varied. This can, for instance, 
explain individual fluctuations spikes, which correspond 
to the crossing of the tunnelling doublet with a third state 
(or an unresolved doublet) 22, FIG. 7] [H FIG. l],|2i 
FIG. 4]. As far as the average of Ae is concerned, we can 
also establish a correspondence between its behaviour as 
a function oil/h and some features of the level dynamics. 
First, a technical point should be mentioned at this stage. 
For r-periodic systems whose classical phase-space is un- 
bounded, it is expected that the quasi-energy spectrum 
becomes a dense set because of the foldings of an infi- 
nite spectrum in the finite Floquet zone ] — 7rfi./T, 7r/i/r]. 
Hence, we need a criterion to select only those states that 
are relevant in the level dynamics. Numerically the trun- 
cation of the Floquet matrices at high |p|'s is not suffi- 
cient, especially for small h, since more and more regular 
and chaotic states are a priori potentially implicated in 
the tunnelling dynamics. To select the levels that are 
actually involved, we will use a simple criterion based on 
the systematic computation of the overlaps am between 
the states li/"™) and a coherent state centered on the is- 
land: cr,„ = K^stabic orbitlV'm)!- Retaining Only thosc 
states whose am is larger than a given threshold (Jfiitor 
will filter the levels that are the best candidates to play 
a role in the tunnelling process, precisely because their 
wave-functions are not negligible in the area where the 
tunnelling doublet wave- functions live. 



Let us fix the classical dynamics at 7 = 0.72 where, 
as we have seen in the previous section, the internal 3/7- 
resonance is suspected to provoke the enhancement of 
average tunneUing for 11 < l/?i < 19 (see FIG. 121. When 



we look at this plateau more carefully, we can identify 
two "bumps" that are visualised in FIG. |3| One cor- 
responds to 11 < 1/h < 13.5 followed by a wider one 
for 13.5 ^ l/fi < 19. These two bumps are in one-to- 
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FIG. 3: (Color online) The upper panels show Aeo as a function of l/ft for 7 = 0.72 and the lower panels display a part of 
the corresponding Floquet spectrum. Only the quasi-energies associated with states whose overlaps are Om > 7 x 10"'^ are 
shown, once calibrated with respect to the tunnelling doublet (unresolved in the lower panels). For relatively large h (on 
the left-hand side), the huge fluctuations of Aeo are in one-to-one correspondance with the crossings of the tunnelling doublet 
by one isolated doublet. On the right-hand side, the upper panel shows a magnification of the plateau in Aeo occurring at 
7 = 0.72 between ft ~ 11 and ft ~ 19 (see also FICl^J. The lower right panel displays a part of the corresponding Floquet 
spectrum. The plateau in Aeo corresponds to a large number of crossings where we can even identify by eye two families of 
levels (encircled by the dashed lines) that give rise to the two bumps in the average of Aeo (red/grey thick continuous line in 
the upper right panel). Beyond the point where the plateau ends, i.e. for ft~^ > 19, no more crossings can be identified, and the 
only other state apart from ej that significantly overlaps with the Husimi function in the island is the bystander doublet e^ 
(see also FlCgJl. 



one correspondence with two bunches of levels with si- 
gnificant overlaps Om, crossing the tunnelling doublet . 
More generally, we have observed in many cases, with 
such a level-dynamics point of view, that there seems to 
be a clear change of regime: for a given CTfiitcr, when 
the average of Ae„ stops decreasing, the states IV'm) 
whose quasi-energies are in the neighbourhood of e,^ 
and whose (T,„ are larger than CTfiitcr become signifi- 
cantly more numerous. For instance, in FIG. 13 where 



In order to be more convincing, this qualitative obser- 
vation can be strengthened by a phase-space analysis. In 
addition to the unresolved tunnelling doublet e^., we can 
clearly see, in the level dynamics in FIG. 01 two other un- 
resolved doublets and for l/h> 20 whose Husimi 
distribution is mainly located near the boundary of the 
stable islands J±. In fact, the doublet can be fol- 
lowed along a quasi-continuous line even for \/h ^ 20 
whereas the other one cannot be identified unambigu- 
ously in that region of h corresponding to the plateau: 



we chose CTfiitcr — 7 X 10~^, or in FIG. 0] with a five times 
more rough filter aeiter — 1-4 x 10~^, just after the end of 
the large plateau at > 20, the crossings with relevant 
levels become suddenly scarce. To put it in a different 
way, the average enhancement of tunnelling appears to 
be the outcome of a collective dynamics involving numer- 
ous states that can be seen in the spectrum through their 
crossing with the tunnelling doublet. 



I 

when 1/h is decreased from about 19, the doublet line 
encounters many avoided crossings and ramifies into the 
bunch of levels we have precisely associated with the 
bump. Therefore, even if their Husimi plots both look 
very similar, one doublet is dramatically involved in the 
tunnelling process while the other remains a bystander. If 
we compute the overlap of these states with the harmonic 
states inside the islands, that is, the eigenstates of the 
harmonic approximation of the Hamiltonian |^ near the 
corresponding stable periodic orbits, we find that the by- 
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FIG. 4: (Color online) The level dynamics for 7 = 0.72 calibrated with respect to the tunnelling doublet (unresolved). 
Compared to FlC^more states are plotted since we take a smaller (Tgitcr — 1.4 x 10~^. When looking at the crossings with 
the tunnelling doublet, there still is a clear transition of regime for below or above 20, which corresponds to the end of the 
plateau in Aeo. The three upper views show the Husimi distribution of the state, \4>t)j whose quasi-energy is et, for n,-^ > 20. 
A comparison with the Poincare surface of section in the neighbourhood of the stable island 1+ is shown. The three lower 
views (which are, actually, very similar to the upper ones at this resolution) correspond to \4>t)j symmetric state whose 
quasi-energy is can be followed for every h~^. Unlike \4>f), the doublet |cjf>^) does not cross and has no influence on the 
tunnelling doublet. (resp. \(f>^)) is localised near the boundary of the two islands and indeed corresponds to the symmetric 

combination of the 7th (resp. 5th) excited regular state in I±. 



stander doublet ef is indeed the 5th excited doublet in X± 
whereas the other ef is the 7th. This can also be checked 
in the Floquet spectrum. Indeed, we expect that the lev- 
els of the local eigenmodes of the island (which is locally 
equivalent to a harmonic oscillator) approximately differ 
from each other by multiples of Hloq where ujo ~ 0.4 de- 
notes the frequency oscillations around the center of the 
island at (p, q) ~ (±1.20, 0). This yields ^ Cq + ^^^o 
modulo a Floquet width 27r/i/r for the levels. The fact 
that it is precisely the quantum number £ = 7 which is 



involved in the emergence of the plateau is not a coinci- 
dence. This £ is exactly the order of the resonance 3/7 
that dominates in X± for 7 = 0.72. 

The quantitative details of how the classical resonances 
may be implemented in order to reproduce the average 
behaviour of Aeo will be described in the next section. 
Even though the approximations that are involved are 
not always under rigorous control, we can see here the 
resonances at work. Suppose that the classical parame- 
ters are such that one classical resonance s/£ dominates 
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the others in the stable islands X± (s being an integer 
and t a strictly positive integer, the order of the reso- 
nance). This precisely means that one torus is actually 
broken into a chain of i sub-islands centered about a sta- 
ble periodic orbit of period 2tt/luq — £t/s. A quantum 
resonance occurs when two quasi-energies em and eg are 
nearly degenerate, i.e. differ by an integer number s' 
of Floquet widths: — eo — s'h^n/r. As mentioned 
above, we also have em — eo — mhujo + s"2tt/t from a 
semiclassical argument (s" being an integer). Hence, we 
immediately see that m must be an integer multiple of 

The analysis given so far, however, shows that the tun- 
nelling enhancement cannot be explained with just one 
crossing of the doublet by e^. This global crossing 
is actually made up of many elementary crossings whose 
contributions cannot be individually distinguished. In ef- 
fect, this game involves many players among which are 
the states that are completely delocalised in the chaotic 
sea. There are also states ("beach" states (41i ,4^] or 
"Janus" states) that are strongly coupled on one side to 
the regular states and on the other side to the chaotic 
ones. It is scarcely surprising that these states have their 
Husimi distribution localised near the borderline between 
the stable islands and the chaotic sea (see FIG. El . 



V. RESONANCE-ASSISTED TUNNELLING 

The above discussion has provided overwhelming evi- 
dence for the relevance of nonlinear resonances in the dy- 
namical tunnelling process. We now focus on the quan- 
titative evaluation of the influence of such resonances, 
which was presented in detail in Refs. ^ IH, [H, |23| . To 
this end, we formally introduce an integrable Hamilto- 
nian that approximately reproduces the dynamics within 
the regular island. For the upper island at p ~ 1, such 
an integrable system can be explicitly obtained by leav- 
ing out the 7+-dependent term in the Hamiltonian 
Performing the time-dependent canonical transformation 
q^q = q — tio the frame that co-propagates with the 
resonant orbit, this integrable Hamiltonian reads 

H,{pA)=^^^- \ cos(g). (3) 

Canonical perturbation theory [3lj | can be applied on the 
basis of Eq. (PJ in order to obtain an improved integrable 
description that is in good agreement with the motion in 
the regular island also at finite values of 7. 

In the dynamics generated by Hq, the regular island 
is embedded in a phase-space domain of bounded ellip- 
tic motion in momentum space. Within this bounded 
domain, action-angle variables (/, 9) can be introduced, 
which respectively correspond to the area enclosed by an 
elliptic invariant orbit as well as to the propagation time 
that elapses along this orbit. In this action- angle vari- 
able representation, we have iJo(p, ?) = Hq{I), and the 
full time-dependent Hamiltonian ^ is formally written 




— TT TT 




FIG. 5: A Janus state: it is strongly coupled to the most 
excited resonant state in the stable island and to the chaotic 
states. In the random matrix model, it is "seen" as chaotic 
from the regular states and almost "regular" from the chaotic 
sea. Here we plot several views of the Husimi distribu- 
tion of a Janus state and compare its localisation with the 
classical phase-space structures (7 = 0.72, = 30, and 

e = -0.007774). 



as H{I,0,t) = Hq{I) + V{I,9,t) where V represents a 
weak perturbation within the island. 

We now assume the presence of a prominent s/£ re- 
sonance within the regular island, where s oscillations 
around the island's center match £ periods of the driv- 
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FIG. 6: Left panel: classical phase space in the vicinity of the 
regular island. The straight solid lines indicate in which order 
the sub-islands are "visited" in the course of time evolution. 
The right panel shows the phase space structure that would 
result from the pendulum-like Hamiltonian Q describing the 
dynamics in the vicinity of the 3/7 resonance. 

ing. This resonance condition is satisfied at £luq{I) — s 
where uq = dHo/dl is the oscillation frequency along the 
bounded orbit with action variable /. The dynamics in 
the vicinity of such a resonance can be approximately 
described by the pendulum-like integrable Hamiltonian 

77eff(/, ^) = ^^—^ + 2Vo cosm (4) 

which is derived from H(I,9,t) using secular perturba- 
tion theory [sj. Here, d = 9 — {s/l)t \s the (slowly vary- 
ing) angle variable that co-rotates with the resonance. 
The resulting phase space structure of H^s is plotted in 
Fig. El 

Following the lines of Ref. j we now evaluate the in- 
fluence of such a resonance in the corresponding quantum 
system by a direct quantization of the effective pendulum 
Hamiltonian in the modified angle variable Apart 
from a phase factor containing the Maslov index, the un- 
perturbed eigenstates of are then given by the plane 
waves {'d\n) = exp(ind). In the co-rotating frame, their 
eigenenergies approximately read 

[ft(n-H/2)-/o]2 

En = ^ (5) 

2mo 

using the fact that /„ = h{n +1/2) are the quantized 
action variables within a regular island of elliptic shape. 
The "i^-dependent term in Eq. Q introduces couplings 
between the states |n) and |n±£) with the coupling ma- 
trix element Vq. These couplings give rise to perturbative 
chains by means of which eigenstates with low and high 
excitations within the bounded domain are connected to 
each other. 

The pendulum Hamiltonian Q can be considered to 
be appropriate for I < Ic where Ic denotes the action 
variable of the outermost invariant elliptic curve of the 
regular island. The regime beyond this "chaos border" 
is characterized by the presence of multiple overlapping 
resonances, which implies that the unperturbed states 
with /„ > Ic can be assumed to be strongly coupled to 
each other by many different matrix elements. Such cou- 
plings would also occur between "bound" and "unbound" 



eigenstates of the integrable Hamiltonian Hq (which are, 
respectively, located within and outside the domain of 
bounded motion that embeds the regular island) as well 
as between states that are located in the vicinity of the 
two different islands at p ~ 1 and p ~ — 1. In this way, 
an efhcient two-step mechanism is introduced by which 
two symmetry-related quasimodes that are localized in 
the upper and lower island, respectively, are coupled to 
each other: the nonlinear resonance connects those quasi- 
modes to the states in the chaotic domain, and the latter 
"see" each other via strong matrix elements of the full 
(Floquet) Hamiltonian of the driven system. 

Fig. El schematically displays the effective Hamiltonian 
matrix that governs this resonance- and chaos-assisted 
tunnelling process between the two "central states" of 
the islands (i.e., given by |n = 0) in the above notation). 
The matrix is restricted to basis states of one particular 
symmetry class — i.e., to states that are "even" or "odd" 
with respect to time-reversal symmetry — and includes, 
for the sake of clarity, only those regular components 
to which the central state is perturbatively connected 
via the s/i resonance. Altogether, kc regular states are 
included in the matrix, where the integer kc is defined 
such that < Ic < Ik^e holds true, i.e., \{kc — 1)£) 

is still located in the island, whereas \kc£) is dissolved in 
the chaotic part. 

We can now prediagonalize the upper left "regular" 
block of the Hamiltonian. This yields, in lowest nonva- 
nishing order in the perturbation strength |Vb| (which is 
much smaller than all relevant energy differences), the 
modified central state as 

This perturbed central now exhibits a nonvanishing ma- 
trix element with one of the states that are contained 
within the chaotic domain: we obtain 

Kff ^ {kcmO) = n E, ^~Em 

k—1 

which can be interpreted as the effective resonance- 
induced matrix element between the central state and 
the chaotic domain. In this way, our effective Hamil- 
tonian matrix can be related to the phenomenological 
matrix models that constitute the starting point of the 
statistical theory of chaos- assisted tunnelling |2lL |2^ l23l | . 

In the simplest possible approach, we assume that the 
"chaos block" is essentially homogeneous and can be well 
modeled by a random matrix from the Gaussian Orthog- 
onal Ensemble (GOE). This assumption does not account 
for the influence of prominent partial barriers to the clas- 
sical transport, which can arise from broken invariant tori 
(so-called "Cantori" j4j|) as well as from unstable perio- 
dic orbits in the chaotic sea, and which would lead to an 
effective division of the chaotic Hamiltonian into several 
sub-blocks that are weakly connected to each other [2l| . 
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FIG. 7: Sketch of the effective Hamihonian matrix that 
describes the coupling between the regular island and the 
chaotic domain for one particular symmetry class (i.e., for 
"even" or "odd" states with respect to time-reversal symme- 
try). The regular part (upper left band) includes only com- 
ponents that are coupled to the island's central state by the 
s/l resonance. In the simplest possible approximation, the 
chaotic part consists of a full sub-block with equally strong 
couplings between all basis states with actions beyond the 
outermost invariant torus of the island. 



Neglecting those partial barriers and performing the ran- 
dom matrix average over the eigenvectors and eigenvalues 
of the chaos block gives rise to the Cauchy distribution 



P{^Eo) = - 



AEo 



TT (AEo)^ + {AEo 



(8) 



for the level splitting between the "even" and "odd" com- 
bination of central states that are associated with the pair 
of regular islands [2^ . This distribution is character- 
ized by the scale 



AEn 



2^ 

NrAr 



(9) 



which contains the most dominant effective matrix ele- 
ment (O between the island's central state and one of 
the chaotic states, the total number Nc of chaotic states, 
as well as the mean level spacing Ac in the chaos block. 
In our case of a periodically driven system, the latter is 
given by = fi/Nc, due to the fact that the chaotic 
states are, in the framework of the Floquet approach, 
uniformly distributed at random in an energy window 
of the size 2'Kh/T = h (we recall that the period of the 
driving equals r = 2ti). As a consequence, we obtain 
A^ = 2^V^^/h. 

The distribution ((SJ is, strictly speaking, valid only 
for AEo VeS and exhibits a cutoff at AEq ~ 2Vcs 
[2^. which ensures that the statistical expectation value 
(AEq) = xP{x)dx does not diverge. However, since 



tunnelling rates and their parametric variations are typi- 
cally studied on a logarithmic scale (see Fig. ^ , we com- 
pute from Eq. JHJ not the mean value (AEq), but rather 
the average of the logarithm of Aii^o- Our "average" level 
splitting (Ai?o)g is therefore defined by the geometric 
mean (Ai?o)g = exp [(ln(A£'o)}] the evaluation of which 
does not involve the above cutoff; we obtain the expres- 
sion 



(AEo), 



AEo 



(10) 



which, notably, is free of any adjustable parameter. 
Hence, up to a trivial prefactor, the mean value of split- 
tings is, in a logarithmic-scale representation, given by 
the square of the coupling matrix element jO between 
the island's central state and the chaos [2£ 
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FIG. 8: Comparison between the numerically calculated 
splittings of the driven pendulum at 7 = 0.72 (black dots) 
and the semiclassical prediction according to Eq. (1101 (red 
solid line). The semiclassical theory, which is based on the 
prominent internal 3/7 resonance within the regular island, 
reproduces quite well the positions and heights of the two 
plateaus that appear in the exact quantum splittings. The 
dashed lines indicate the size of the logarithmic standard de- 
viation according to Eq. IIH . The green dot-dashed line dis- 
plays the prediction that would be obtained from Eq. H12|l 
with Q = 1. 

Figure |S1 shows the comparison with the exact quan- 
tum splittings at 7 = 0.72, calculated by the numerical 
diagonalization of the Floquet matrix. The semiclassi- 
cal prediction H10|l (solid line in Fig. IHJ was evaluated 
on the basis of the prominent 3/7 resonance, for which 
the relevant parameters /q, mo and Vq that enter into 
the pendulum Hamiltonian (@J) were entirely determined 
from classical dynamics of the system: As in Ref. ,,2 6.] , we 
compute for this purpose the trace of the monodromy ma- 
trix associated with a stable or unstable periodic point 
of the 3/7 resonance, as well as the phase space areas 
that are enclosed by the inner and outer separatrices of 
the resonance. Indeed, those quantities remain invariant 
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under the canonical transformation to the action-angle 
variables {!,"&), which means that the latter need not be 
explicitly evaluated in order to obtain the effective cou- 
pling matrix element VcS- 

We see that the semiclassical theory reproduces quite 
well the two plateaus that arise in the quantum splittings. 
The drop in the semiclassical splittings at 1/fi ~ 24 oc- 
curs due to the fact that the island supports more than 
seven locally quantized eigenstates beyond this critical 
value of 1/h; hence, two perturbative steps instead of 
one are required in order to connect the central state of 
the island to the chaotic domain. This drop is consider- 
ably softened in the quantum splittings, but can still be 
identified, which confirms the relevance of the resonance- 
assisted coupling mechanism in this tunnelling process. 
Since nonlinear resonances represent a general feature of 
nonintegrable Hamiltonian systems, we expect that the 
appearance of such step-like sequences of plateaus is a 
generic phenomenon, which is not restricted to one par- 
ticular system, but arises in various chaos-assisted tun- 
nelling processes. This expectation is indeed confirmed 
by previous studies on the kicked Harper on the 

decay of nondispersive wave packets in driven hydrogen 
|45|. as well as on the dynamical tunnelling process of 
cold atoms 0, where significant plateau structures 
were encountered in numerically calculated tunnelling 
rates. 

In addition to the mean value for the splittings, the 
probability distribution (jSJ can also be used in order to 
calculate the expectation value for their logarithmic vari- 
ance characterizing the average size of fluctuations on a 
logarithmic scale; we obtain 



(11) 



This result is universal in the sense that it does not de- 
pend on system-specific parameters nor on h. The size of 
the corresponding standard deviation is indicated by the 
dashed lines in Fig. |H1 which are generated by multiply- 
ing AEq with e^^^^. Clearly, the "window" defined by 
those dashed lines is exceeded near local avoided cross- 
ings with chaotic states, where a large enhancement as 
well as a complete suppression of the splittings may be 
induced Apart from those exceptional events, how- 
ever, the scale of the average fluctuations of the splittings 
is well described by Eq. fTT|l . The green dot-dashed hue 
in Fig. |S1 displays the prediction for the splittings that 
would be obtained from Eq. jSJ according to Ref. [2^ . 
Here, the asymptotic expression 



Ae 



-2(l-ln 2)N 



(12) 



of Eq. |(2Jl, vahd for N = A/{2Trh) > 1, was evaluated 
with the numerically calculated size A of the island. Be- 
ing an entirely classical (fi-independent) quantity, the 
system-dependent prefactor O was set equal to unity. 
We clearly note a substantial disagreement between this 
semiclassical estimation and the quantum splittings. 



The theory of resonance-assisted tunnelling does not 
reproduce the splittings in the deep quantum regime of 
larges values of h, where a different mechanism, possibly 
in the spirit of Ref. [23|. might induce the transition to 
the chaos. Moreover, the critical value of 1 /h at which the 
drop of the splittings from the first to the second plateau 
occurs is significantly overestimated by our semiclassical 
approach. This could be due to the presence of partial 
barriers in the chaotic phase space domain, such as "Can- 
tori" [4^ . which are known to inhibit the quantum flux 
at not too small values for h 0, ■ 1^ such a Cantorus 
is manifested in the vicinity of the regular island, the 
effective "quantum" size of the island could be consid- 
erably enhanced as compared to A, which would reduce 
the value of 1/h at which exactly £ quantum states are 
localized around the island. This observation is indeed in 
accordance with the manifestation of Janus states [ZHIi^ 
in the spectral analysis (see Section Hvfl . 

Significant deviations of the splittings from the semi- 
classical prediction IjlOl) are also to be expected in the 
deep semiclassical regime where a multitude of steps 
{k 3> 1) would be needed to connect the central state 
to the chaotic domain according to the expression 10. 
In this regime, the coupling via the s/i resonance — 
which also represents a dynamical tunnelling process as 
was pointed out in Ref. — can again be assisted 
by the presence of another nonlinear s' /£' resonance, as 
long as this is permitted by the associated selection rule 
{n t-^ n + i'). Such a s' /£' resonance would generally 
exhibit a lower effective coupling strength Vq and is typi- 
cally of higher order than the s/£ resonance (i.e., £' and s' 
are typically larger than £ and s, respectively). The rel- 
evance of this multi-resonance coupling mechanism was 
demonstrated in the near-integrable kicked Harper model 
where the semiclassical tunnelling process involves a se- 
quence of three nonlinear resonances 0) HE 113 ■ 




FIG. 9: Classical phase space at 7 = 0.67 in the vicinity of the 
regular island I+. The two chains of sub-islands correspond 
to the 3/7 as well as to the 5/11 resonance (inner and outer 
chain, respectively). 

In the case of the driven pendulum at 7 = 0.72, such 
multi-resonance processes cannot be observed within the 
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range of values for 1/h at which precise quantum cal- 
culations of the splittings can be performed. This is 
different, however, for 7 = 0.67. At this value of the 
coupling parameter, the regular island T± exhibits two 
nonlinear resonances of almost equal importance: a 3/7 
resonance, located closer to the center of the island than 
at 7 = 0.72, and a 5/11 resonance, located close to the 
chaos border (see Fig. This implies that a two-step 
tunnelling process involving both resonances can be en- 
countered at finite values for \/h. Assuming that each 
resonance contributes with one single perturbative step 
(which is the case at the values of 1/S considered here) , 
the corresponding matrix element that connects the cen- 
tral state to the chaotic domain would read 



Vr 



(3/7) 



(3/7) 



E. 



(3/7) 



(5/11) 



(13) 



has the larger matrix element and should there- 

fore dominate compared to the 3/7 resonance) and for 
l/h > 28.8 obtained through the two-step process that 
is described by Eq. (|13|l . As in the case of 7 = 0.72, 
the overall agreement between the semiclassical and the 
quantum splittings is quite good, with significant de- 
viations arising only in the deep quantum regime at 
l/h < 10 as well as in the vicinity of the cross-over be- 
tween the single- and the two-step process, which is ar- 
tificially sharp in the semiclassical calculation. Striking 
evidence for the validity of the two-step process is the ap- 
pearance of a pronounced peak in the quantum splittings 
at \/h ~ 43, which arises due to a vanishing denomina- 
tor in Eq. (|13|l . i.e., due to the fact the central state and 
the 7th excited state become near-degenerate in the co- 
rotating frame. The position and height of this peak are 
fairly well reproduced by the semiclassical theory. 



Here, Vq represents the pendulum coupling strength 
of the effective Hamiltonian and En^^'^ denote the en- 
ergies (ISJ in the co-rotating frame that is defined with re- 
spect to the sj I resonance. This expression requires that 
the 7th excited state, i.e. the first state to which the cen- 
tral state of the island is coupled via the 3/7 resonance, 
is located in between the two resonances, which implies 
that the action variable /y of this state ought to be lower 



than the action variable /, 



(5/11) 



of the 5/11 resonance. 



This condition turns out to be valid for 1/fi, > 28.8 



10 



10"' 



10" 



10" 




FIG. 10: Comparison between the numerically calculated 
splittings of the driven pendulum at 7 = 0.67 (black dots) and 
the semiclassical prediction (red solid line). The latter was 
evaluated by a single-resonance process via the 5/11 resonance 
for 1/ft < 28.8, and by a two- resonance process according to 
Eg. 1131 involving also the 3/7 resonance for 1/ft > 28.8. As 
in Fig. |H| the dashed lines indicate the size of the logarithmic 
standard deviation according to Eq. lllll . and the green dot- 
dashed line displays the prediction obtained from Eq. 1121 . 

Fig. ^1 displays the quantum splittings at 7 = 0.67. 
The red solid line denotes the semiclassical prediction 
(|10|l for the splittings, which is for l/h < 28.8 calculated 
by the single-step process via the 5/11 resonance (which 



VI. CONCLUSION 

In summary, we provided clear evidence of the signifi- 
cance of nonlinear resonances in the dynamical tunnelling 
process within the driven pendulum. Indeed, the sig- 
nature of the 3/7 resonance that is dominantly mani- 
fested in the regular island at 7 = 0.72 is identified in 
various ways: (a) The characteristic plateau structure 
in the level splittings sensitively depends on whether or 
not the resonance is well developed within the regular 
island, (b) Floquet states that exhibit an appreciable 
overlap with the central Husimi wavefunction of the is- 
land energetically correspond to the 7th excited eigen- 
state within the island, and (c) a semiclassical expres- 
sion for the tunnelling-induced splittings that is based on 
the resonance-assisted coupling scheme is in good agree- 
ment with the exact quantum data. The validity of the 
resonance-assisted tunnelling mechanism is furthermore 
confirmed by a striking peak in the quantum splittings at 
7 = 0.67, which arises due to a two-step process involv- 
ing two different nonlinear resonances within the regular 
island. Both in the spectral analysis and in the semiclas- 
sical comparison, we identify traces of Janus states that 
are located in the chaotic vicinity of the islands. This in- 
dicates that partial barriers in the chaotic domain could 
still be relevant in our system at the values of h under 
consideration. 

The present approach is presumably not suited for a 
full-blown semiclassical theory of mixed regular-chaotic 
tunnelling in terms of complexified orbits, for which the 
ansatz of Refs. 

El in LfiJ provides a more convenient 
framework. Our findings, however, provides essential in- 
gredients to the interpretation of such semiclassical the- 
ories, in the sense that "direct" and "resonance-assisted" 
processes ought to be somehow represented in relevant 
combinations of such complex orbits. Moreover, specific 
quantitative predictions for tunnelling rates, using only 
easily accessible quantities of the classical dynamics of 
the system, can be made on the basis of our approach. 
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provided h is small enough for the most dominant reso- 
nance to become relevant (which is roughly the case if ^/2 
states fit into the island). We therefore believe that the 
principle of resonance-assisted couplings will represent 
the relevant paradigm also in the context of more com- 
plex dynamical tunnelling processes, e.g., within systems 
that have two or more degrees of freedom. 

Acknowledgments: It is a pleasure to thank Denis 



Ullmo for sharing with us his extensive experience on 
tunnelling. A.M. acknowledges the generous hospital- 
ity of Dominique Delande at Laboratoire Kastler-Brossel 
and his relevant remarks after reading the first proof of 
this manuscript. P.S. acknowledges support from the 
Bayerisch-Franzosisches Hochschulzentrum (BFHZ) and 
from the Deutsche Forschungsgemeinschaft (DFG). 



[1] W. A. Lin and L. E. Ballentine, Phys. Rev. Lett. 65, 
2927 (1990). 

[2] F. Grossmann, T. Dittrich, P. Jung, and P. Hanggi, Phys. 

Rev. Lett. 67, 516 (1991). 
[3] A. Mouchet, C. Miniatura, R. Kaiser, B. Gremaud, and 

D. Delande, Phys. Rev. E 64, 016221 (2001). 
[4] W. K. Hensinger, H. HafFner, A. Browaeys, N. R. Hecken- 

berg, K. Helmerson, C. McKenzie, G. J Milburn, W. D 

Phillips, S. L. Rolston, H. Rubinsztein-Dunlop and B. 

Upcroft, Nature 412, 52 (2001). 
[5] D. A. Steck, H. O. Windell, and M. G. Raizen, Science 

293, 274 (2001). 
[6] W. K. Hensinger, A. Mouchet, P. S. Julienne, D. De- 
lande, N. R. Heckenberg and H. Rubinsztein-Dunlop, 

Phys. Rev. A (3) 70, 013408 (2004). 
[7] C. Dembowski, H.-D. Graf, A. Heine, R. Hoflerbert, H. 

Rehfeld and A. Richter, Phys. Rev. Lett. 84, 867 (2000). 
[8] R. Hoflerbert, H. Alt, C. Dembowski, H.-D. Graf, H. L. 

Harney, A. Heine, H. Rehfeld and A. Richter, Phys. Rev. 

E 71, 046201 (2005). 
[9] L. D. Landau and E. M. Lifshitz, Quantum Mechanics 

(non relativistic theory), Vol. 3 of Course of Theoretical 

Physics (Pergamon Press, Oxford, 1958). 
[10] R. Balian and C. Bloch, Ann. Physics 84, 559 (1974). 
[11] M. Kus, F. Haake, and D. Delande, Phys. Rev. Lett. 71, 

2167 (1993). 

[12] P. Leboeuf and A. Mouchet, Phys. Rev. Lett. 73, 1360 
(1994). 

[13] R. Scharf and B. Sundaram, Phys. Rev. E 49, R4767 

(1994) . 

[14] S. C. Creagh and N. D. Whelan, Phys. Rev. Lett. 77, 
4975 (1996). 

[15] For multidimensional non-separable but integrable sys- 
tems, these tori are precisely the relevant classical struc- 
tures that allow to generalise jwkb Id-tunnelling j3S| . 

[16] A. Shudo and K. S. Ikeda, Phys. Rev. Lett. 74, 682 

(1995) . 

[17] A. Shudo and K. S. Ikeda, Physica D 115, 234 (1998). 
[18] A. Shudo, Y. Ishii, and K. S. Ikeda, J. Phys. A 35, L225 
(2002). 

[19] O. Bohigas, M.-J. Giannoni, and C. Schmit, Phys. Rev. 

Lett. 52, 1 (1983). 
[20] O. Bohigas, D. Boose, R. E. de Carvalho, and V. Mar- 

vuUe, Nuclear Phys. A 560, 197 (1993). 
[21] O. Bohigas, S. Tomsovic, and D. Ullmo, Phys. Rep. 223, 

43 (1993). 

[22] S. Tomsovic and D. Ullmo, Phys. Rev. E 50, 145 (1994). 
[23] F. Leyvraz and D. Ullmo, J. Phys. A 29, 2529 (1996). 
[24] O. Brodier, P. Schlagheck, and D. Ullmo, Phys. Rev. 

Lett. 87, 064101 (2001). 
[25] O. Brodier, P. Schlagheck, and D. Ullmo, Ann. Physics 



300, 88 (2002). 

[26] C. Ehschka and P. Schlagheck, Phys. Rev. Lett. 94, 
014101 (2005). 

[27] P. Schlagheck, C. Eltschka, and D. Ullmo, in Progress 
in Ultrafast Intense Laser Science I, edited by K. Ya- 
manouchi, S. L. Chin, P. Agostini, and G. Ferrante 
(Springer, Berlin, 2006), p. 107. 

[28] V. A. Podolskiy and E. E. Narimanov, Phys. Rev. Lett. 
91, 263601 (2003). 

[29] A. Mouchet and D. Delande, Phys. Rev. E 67, 046216 
(2003). 

[30] B. Chirikov, Phys. Rep. 52, 263 (1979). 

[31] A. J. Lichtenberg and M. A. Lieberman, Regular and 
Stochastic Motion, Vol. 38 of Applied Mathematical Sci- 
ences (Springer- Verlag, New York, 1983) . 

[32] The only difference is that a third resonance term must 
be added to Q in order to fulfil some experimental con- 
straints on the lasers configuration. How this can be done 
is explained in Refs. 0, but the present paper, de- 
voted to theoretical considerations only, will not intro- 
duce this irrelevant third term. 

[33] M. J. Davis and E. J. Heller, J. Phys. Chem. 85, 307 
(1981). 

[34] M. Latka, P. Grigolini, and B. J. West, Phys. Rev. E 50, 
596 (1994). 

[35] V. Averbukh, N. Moiseyev, B. Mirbach, and H. J. Korsh, 
Z. Phys. D 35, 247 (1995). 

[36] L. Bonci, A. Farusi, P. Grigolini, and R. Roncaglia, Phys. 
Rev. E 58, 5689 (1998). 

[37] For the same reason we will not consider the different 
possible boundary conditions offered by Bloch's theorem. 
The wave-functions in the q representation will always 
be taken to be strictly 2tt spatially periodic, i.e. with- 
out acquiring any phase-factor when winding around the 
cylinder. 

[38] S. C. Creagh, J. Phys. A 27, 4969 (1994). 

[39] R. Egydio de Carvalho and A. P. Mijolaro, Phys. Rev. E 

70, 056212 (2004). 
[40] R. Roncagha, L. Bonci, F. M Izrailev, B. J. West, and P. 

Grigolini, Phys. Rev. Lett. 73, 802 (1994). 
[41] E. Doron and S. D. Frischat, Phys. Rev. Lett. 75, 3661 

(1995). 

[42] S. D. Frischat and E. Doron, Phys. Rev. E 57, 1421 
(1998). 

[43] A. M. Ozorio de Almeida, J. Phys. Chem. 88, 6139 
(1984). 

[44] R. S. MacKay, J. D. Meiss, and I. C. Percival, Phys. Rev. 

Lett. 52, 697 (1984). 
[45] K. Hornberger and A. Buchleitner, Europhys. Lett. 41, 

383 (1998). 

[46] For this latter case, we have indications that the promi- 



13 



nent plateau structure in the level splittings at 7 = 0.18 
(see Fig. 7 in Ref. 0) is induced by a nonlinear 1/4 reso- 
nance that is already contained within the chaotic part 
of the phase space. Since our numerical approach to cal- 
culate the coupling matrix element Vq is perturbative by 
nature and relies on a well-preserved separatrix structure 
of the resonance, we are, at present, unable to make a 



quantitative semiclassical reproduction of this particular 
plateau. 

[47] T. Geisel, G. Radons, and J. Rubner, Phys. Rev. Lett. 

57, 2883 (1986). 
[48] N. T. Maitra and E. J. Heller, Phys. Rev. E 61, 3620 

(2000). 



